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Abstract 

An introduction is given to the many-body coherent potential approximation (CPA) for 
the double-exchange (DE) model and the Holstein-DE model, the latter including coupling 
of the electrons to local phonons as well as to the local spins. It is shown how the method 
can treat the local spins and phonons quantum-mechanically and how it is equivalent to 
dynamical mean field theory in the classical limit. In the Holstein-DE model a full discussion 
is given of the cross-over from weak electron-phonon coupling through intermediate coupling, 
where small-polaron bands begin to appear, to strong coupling where some results similar to 
those of standard small-polaron theory are recovered. The theory is applied to ferromagnetic 
manganites with a full discussion of magnetic, transport, and spectroscopic data. It is found 
that many manganites are in the critical regime on the verge of small-polaron formation, 
which explains their sensitivity to parameters such as applied magnetic field and pressure. 
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1 Introduction 



Magnetic materials exhibiting metallic behaviour can often be considered as systems of local mo- 
ments coupled to electrons in a conduction band by local exchange interactions. The Hamiltonian 
for such a system is 

h de - £ t ijC lc ja -J^Si-n-h^ {Si + Of) , (1) 

ija i i 

where c\ a creates an electron of spin a on lattice site i, Si is a local spin operator and er; = 
(of , of , of) is a conduction electron spin operator defined by 

cxf = of + iof = 4 T c a , or = of - zo? = cJ iCiT j of = — (n lT - n a ) (2) 

with rijcr = c\ a Ci a . The three terms of equation (|l|) describe hopping of the conduction electrons, 
exchange coupling between local and itinerant spins and coupling to an external magnetic field. If 
the local exchange coupling arises from hybridization between the localized and itinerant electrons, 
as in anomalous rare earth systems exhibiting heavy fermion behaviour, the exchange parameter 
J is negative. The Hamiltonian (|lj) is then often called the Kondo lattice model in view of its 
connection with the Kondo impurity model which has a local spin on one site only |Q. When 
Hund's rule coupling is dominant J > and the system is sometimes called a ferromagnetic 
Kondo lattice. This is misleading since for J > there is no connection with the Kondo effect. 

For J > it is useful to distinguish two distinct physical regimes, depending on the magnitude 
of J compared with the width 2W of the conduction band. If J <C W, as in a normal rare earth 
metal, the exchange coupling can be treated as a perturbation which gives rise to the Ruderman- 
Kittcl-Kasuya-Yosida (RKKY) interaction between local moments. In most rare earth metals this 
interaction, which oscillates in space, leads to oscillatory or spiral configurations of the localized 
f electron moments. The uniform ferromagnet Gd is an exception. In this weak coupling regime 
the Hamiltonian (|l|) is usually referred to as the s — f or s — d model. 

If J 3> W the exchange coupling can no longer be treated as a perturbation. A conduction 
electron can only hop onto a site with its spin parallel to the local moment at that site. Furthermore 
if the number of conduction electrons per atom n < 1 double occupation of a site is strongly 
suppressed. A single electron at a site, with its spin parallel to the local spin S, enjoys an 
exchange energy —JS/2 which is lost if a second electron hops on. The system is therefore a 
strongly correlated electron system, just like the Hubbard model in the regime of strong on-site 
Coulomb interaction U , and for n — 1 the system is a Mott insulator. In much of the theoretical 
work on the present model the local spins are treated as classical vectors, corresponding to S — > oo. 
Since for J 3> W the itinerant spin must always be parallel to the local spin on each site, the 
effective hopping integral for hopping between sites i and j becomes t { - cos (%/2j, where 9^ is 

the angle between the classical spins Sj. The cosine factor arises from the scalar product 
of two spin 1/2 eigenstates with different axes of quantization. The resultant band narrowing 
in the paramagnetic state favours ferromagnetism in order to lower the kinetic energy. This 
mechanism for ferromagnetism was first introduced by Zener Q and developed by others || |[ ||] . 
Since it involves strong exchange coupling on two adjacent atoms it is known as double-exchange. 
Consequently the Hamiltonian (|l|) in the strong-coupling regime J » W is called the double- 
exchange (DE) model. This paper is concerned with the DE model, with quantum and classical 
local spins, and with its extension to include coupling of the electrons to local phonons. We call 
this extended model the Holstein-DE model. 

The paper provides a simple introduction to the many-body coherent potential approximation 
(CPA) described in detail in our previous papers j^, ^, ||, ^| and reviews the results. When spins and 
phonons are treated classically the CPA is equivalent to the dynamical mean field theory (DMFT) 
approach of Furukawa jl(| [ll| and Millis et al. |12] |. A fully quantum mechanical treatment within 
DMFT is a difficult numerical task which has not yet been done. The many-body CPA may be 
regarded as a useful analytic approximation to DMFT which becomes exact in the classical limit. 
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The main application we have in mind is to the manganites La lx D x Mn0 3 , where D is a 
divalent ion such as Ca or Sr. Recently there has been renewed interest in these systems due to the 
discovery of colossal magnetoresistance (CMR), the name given to a large reduction in resistivity 
in an applied magnetic field [ p^[ . For x = 0, LaMn0 3 is an antiferromagnetic insulator, but on 
doping with 0.2 < x < 0.5, the system becomes a ferromagnetic metal. The CMR phenomenon is 
observed in the vicinity of the Curie temperature T c . In applying equation ([!]) to these systems 
the local spins S i are of magnitude S — 3/2, corresponding to three localized Mn d electrons 
of t 2g symmetry, and the band is derived from Mn d states of e g symmetry. The band contains 
n = 1 — x itinerant electrons per atom. Since we are interested in doped ferromagnetic systems, 
the antiferromagnetic interaction between neighbouring local spins plays no important role and 
is neglected. A more serious simplification is the use of a single s band, instead of two d bands 
based on orbitals of e g symmetry. We discuss this in section |?]. 

In section |^ we introduce the many-body CPA by means of the Hubbard model and in section 
^ we develop it for the more complicated case of the DE model with quantum local spins. In 
section [| we discuss the electrical resistivity of the paramagnetic state in the DE model and 
show, as first stressed by Millis et al. |l4[, that the DE model is unable to describe the physics 
of the manganites completely. In section |5j we briefly discuss the variation of T c as a function of 
band-filling n. The development of the many-body CPA for the Holstein-DE model is discussed 
in section | and a comparison with standard small-polaron theory is made in the limit of strong 
electron-phonon coupling. For intermediate coupling we discuss the crossover, with increasing 
temperature, from polaronic behaviour to a situation where the phonons behave classically, the 
case considered by Millis et al. In section [?| we consider the application of the Holstein-DE 
model to the manganites. Comparison with experimental data leads to the conclusion that typical 
manganites lie in the critical intermediate coupling regime which is not fully described by previous 
theories. We summarise briefly in section H. 



To introduce the many-body CPA we consider the Hubbard model, which is a simpler model for 
strongly correlated electrons than the DE model with quantum spins. The Hamiltonian for this 



and Hubbard |15| set out to find an approximation to the one-electron retarded Green function 
G ka (e) which is exact in the atomic limit — 0. He used the equation of motion method and 
the idea of the alloy analogy described below. It turns out that Hubbard's approach, without 
the minor "resonance broadening correction" , is equivalent to the CPA which was developed later 
p6| . The CPA derivation of Hubbard's result is much simpler than the original equation of motion 
method. However we had to resort to an extension of the original method to derive the many-body 
CPA for the DE model, with and without phonons, in the general form needed to discuss magnetic 
properties. In this paper we restrict the derivation to the paramagnetic state in zero magnetic 
field, although we summarize some more general results. 

The alloy analogy consists in considering the j-spin electrons, say, to move in the potential 
due to static |-spin electrons, frozen in a random configuration which must be averaged over. 
Thus a one-electron Hamiltonian for |-spin is obtained from by taking the last term to be 
U J2i n it( n ii) where {n^) takes the value 1 with probability and with probability 1 — n^. 
Here n a is the number of er-spin electrons per atom. It is important to note that the alloy analogy 
is quite distinct from the Hartree-Fock approximation in which (n-) = n a for all i. In the alloy 
analogy a a-spin electron moves in a random potential given by U on n~ sites and on 1 — n- 
sites. In the CPA the random potential is replaced by a uniform, but energy-dependent, effective 
potential S_(e) for an "effective medium". This effective potential, in general complex, is called 
the coherent potential and is in fact the electron self-energy. The procedure for determining 



2 CPA for the Hubbard model 



model is 




(3) 
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is to insist on a zero average t-matrix for scattering by a central atom, with potential U or 0, set 
in the effective medium. Equivalently the average of the site-diagonal element G (e) of the Green 
function, for each type of central atom, is put equal to the site-diagonal element of the Green 
function for the effective medium. Thus 

G - = n * i-(u G -\ )Ga + ( J - »*> T+rfe w 

v a > a a a 

and 

where the bare band Green function is given by 

G (e) = / de'^ . (6) 

Here e k — J2j tij exp[ifc • (R i — R.)] is the band energy, where R^ is the position of site i, N (e) is 
the corresponding density of states per atom and N is the number of lattice sites. Equations (Q) 
and (||) are to be solved self-consistently for S CT (e) and hence for the local Green function G CT (e). 
Equation (|j) may be written as 

G Ug \ 1 ~ H ° (7) 

This may be compared with the exact Green function for the atomic limit (t- — 0) which is given 
by@ ' 

^(e) = ^ + ^, (8) 

where in this retarded Green function e has a small positive imaginary part. Hence 

G tr (e)=G? (S. + G; 1 ) . (9) 

Clearly this CPA equation is exact in the atomic limit, when N (e) = 6(e) and it follows from 
equation (|J) that S CT + G^ 1 = e. Solution of the CPA equation becomes simple if the density of 
states N Q (e) is taken to be of the elliptic form 

W = ^(^ 2 - 2 ) 1/2 do) 



where 2W is the bandwidth. Then from equation (Q) 
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G Q (e) = ^[e-(e^W^ 1/2 ) . (11) 

Introducing this expression for G a in equation (|^), and solving for e — S CT (e), we find 

K + G- 1 = e - W 2 GJA (12) 

Hence equations (||) and (||) give an algebraic equation for G a . 

Solving this type of equation for G a , Hubbard Q calculated the density of states N a (t) = 
— 7r ImG CT (e), considering particularly the paramagnetic state = n^ = n/2 and concentrating 
on the half- filled band case n = 1. He showed that for U/W greater than a critical value, equal to 
1 in the present approximation, a gap opens in N(e) at the Fermi level so that the system becomes 
an insulator as envisaged by Mott. For U ^S> W the density of states consists of two peaks centred 
on e = and e = U, these being broadened versions of the ^-functions at these energies in the 
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atomic limit. Furthermore, for general band-filling n, the spectral weights in the two peaks are 
the same as in the atomic limit. 

The CPA for the Hubbard model has some serious defects. There are no self-consistent solutions 
with magnetic order. Furthermore in the paramagnetic metallic state, for n < 1 or for n = 1 with 
U/W less than the critical value, the system is never a Fermi liquid. There is never a sharp Fermi 
surface at T = with a Migdal discontinuity in the Bloch state occupation number, as pointed 
out by Edwards and Hewson |T^] . This is due to the absence of states with infinite lifetime at the 
Fermi level, since within the alloy analogy all states are scattered by disorder. A modification of 
the CPA to remedy this defect, retaining the analytic simplicity of the method, had some limited 
success Jl8|, [n], |(| . However the most satisfactory approach is DMFT which involves numerical 
solution of an associated self-consistent impurity problem [^l], DMFT may be regarded as 
the best local approximation, in which the self-energy is a function of energy only, and is exact in 
infinite dimension. 

The many-body CPA is considerably more satisfactory for the DE model than it is for the 
Hubbard model, as discussed in the next session. There is one limit, the case of classical spins 
(S = oo), in which the CPA is identical to DMFT. This is because classical spins are static and 
an alloy analogy of frozen disordered spins is completely justified. DMFT for the DE model has 
only been implemented for classical spins [ll[] and the many-body CPA discussed in the next 
section provides an approximate analytic extension of DMFT to quantum spins. The system 
orders ferromagnetically below a Curie temperature T c , as it should, and the disordered spin 
state above T c should be well described. However, the accuracy of the ground state at T = 
for finite S is unclear. The saturated state with all itinerant and local spins completely aligned, 
which is the ground state for S = oo (we are always considering large J in the DE model) , is never 
a self-consistent CPA solution for finite S Q. Actually the parameter range of stability of the 
saturated ground state is unknown. It has been shown rigorously that for J — oo it is unstable for 
S = 1/2 and 0.12 < n < 0.45, with a simple cubic nearest-neighbour tight-binding band (2j|. If 
the true ground state is not saturated it seems unlikely to be a uniform (spatially homogeneous) 
ferromagnet, with partially ordered local and itinerant spins, as in the uniform CPA ground state 
for finite S. Such a state would probably not be a Fermi liquid, just as in CPA, unless the electrons 
making up the spin S became partially delocalised with spectral weight at the Fermi level. Such 
speculation goes beyond the DE model. 



3 The many-body CPA for the DE model 

In an earlier paper || we developed the many-body CPA for the DE model using an extension of 
Hubbard's equation of motion method. Hubbard's "scattering correction" becomes more compli- 
cated owing to the form of the interaction term in the DE model whereby electrons can flip their 
spin via exchange of angular momentum with the local spins. This dynamical effect couples the 
equations for G+ and G, and was first treated by Kubo in a one-electron dynamical CPA. 
The main feature of our many-body CPA is that we recover Kubo's one-electron CPA as n — > 
and the correct atomic limit for general band-filling n as — > 0. In a second paper Q we showed 
the equivalence to DMFT in the limit S — > oo, J — > oo. The full equation of motion derivation 
of the many-body CPA is required to obtain general results in the presence of a magnetic field 
and/or magnetic order jjj. However it turns out that in the zero-field paramagnetic state we can 
deduce the CPA equation from the atomic limit Green function G^ h and equation (||), just as in 
the Hubbard model. We shall therefore not repeat the full derivation in this paper although we 
shall discuss results on magnetic properties in section |L 

The atomic limit Green function, G^ L , say, is easily obtained by the equation of motion method 
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using the Hamiltonian ([!]) with t i 



The result for zero field (h = 0) is 



1 



25 + 1 



(S' + S' z )n i -S , ^ct+ 
e + J(5 + l)/2 



(S-S z )(l-n 



5"^ 



+ 



(S+l + S z )(l-n,) +S-a+ 



e- J(5 + l)/2 
(S'+l-S' z )n i + S , -cr + 



J5/2 



J5/2 



(13) 



and for h ^ one merely has to replace e by e + h/2. The angle brackets (...) represent thermal 
averages and all operators within them correspond to the same site i, this suffix thus being omitted. 
This expression, with four poles, is considerably more complicated than the two-pole Hubbard 
model expression of equation (ph. The poles at e = +J5/2, ±J(5 + l)/2 correspond to energies 
to add or remove an electron from the atom, that is to transitions between singly-occupied states 
and either unoccupied or doubly-occupied states. The singly-occupied states have total spin 5 + ^ 
or 5 — h with energies — JS/2 and J(5 + l)/2 respectively; the unoccupied and doubly-occupied 
states have zero energy. 

In the zero-field paramagnetic case it turns out that the CPA equation for G(e) with the 
redundant suffix a omitted, is given by equation (|J) as in the Hubbard model. Thus, taking the 
band to have the elliptic form (0), the CPA equation for G is 



with G AL given by 



G^(e) 



1 



G(e) = G (e - W 2 G/4) 



iS/2-(S-a) 5(1 - n/2) - (S ■ a) 



(14) 



25+1 



e + J(5 + l)/2 



J(5+l)/2 



(5 + 1)(1 - n/2) + {S-a) n(S + l)/2 + (S ■ a) 



JS/2 



JS/2 



(15) 



The spin symmetry of the paramagnetic state has been used to simplify the expectations in the 
previous form of G AL , equation (|l3|). It is easy to show that {S ■ er) — > nS/2 as J — > oo and (S ■ a) 
will be very near this limit as long as JS > 2W. We make this approximation in calculating G(e), 
and hence the density of states, N(e) — — 7r _1 ImG(e) from equations ( |l4| ) and (p^|). The results 
are shown in figure for S = 3/2 and J = AW for various n. Clearly, from equation (|l5|), the 
approximation to (S ■ <x) has the effect of removing the weak band centred on f = —J(S + l)/2 
but it does not affect the total weight or the distribution of weight between the two lower and 
two upper bands. It may be seen that as n increases from the band near e = J(S + l)/2 is 
reduced in weight and a new band appears near JS/2, until at n = 1 no weight remains in the 
band near J(S + l)/2. The weight in the band near - JS/2 is (S + 1 - n/2)/{2S + 1) per spin 
so if JS is sufficiently large to separate the bands (JS > 2W) this band will just be filled at 
n = 1 producing a Mott insulator as expected. This redistribution of weight between bands as 
they fill with electrons is characteristic of the many-body CPA p4j and was missing from Kubo's 
one-electron CPA which was restricted to n — 0. 

In the strong-coupling limit J — ► oo, which is taken with a shift of energy origin e — > e— JS/2, 
equation (Ea) simplifies to 



G AL (e) 



L (5+l-n/2)/(25+l). 



Equation (p 



then becomes a quadratic equation for G with solution 

, 2 



G(e) 



D 2 



D 2 



(16) 



(17) 



where a 2 = (S + 1 — n/2)/(2S + 1) and D = aW. By comparing with equations ( |l0| ) and ( |TT| ) 
we see that the density of states is a single elliptical band of weight a 2 and bandwidth 2aW. As 
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Figure 1: The density of states in the paramagnetic state of the double-exchange model for 
S = 3/2, J = AW, and n = 0, 0.25, 0.75 and 1. Energy units of W are used. 
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S — > oo the band- narrowing factor a — > 1/V2 = 0.707, which is close to the classical result of 2/3, 
obtained by averaging cos (9/2) over the solid angle. 

In the classical spin limit S — > oo we rescale J, replacing it by J/S, and the CPA equation 
becomes 



G = 



1 



1 



1 



T. + G- 1 + J/2 S + G- 1 - J/2 



(18) 



Here we have used the general equation (g), valid for arbitrary band-shape, rather than equa- 
tion (|lj). Equation (|l^) is precisely the equation obtained by Furukawa |S within DMFT. 



4 Resistivity in the paramagnetic state of the DE model 

The Kubo formula for the conductivity a involves the two-particle current-current response func- 
tion. However in the local approximation of CPA or DMFT there is no vertex connection and a 
may be expressed in terms of the one-particle spectral function 

A fe (e) = -7r- 1 ImG fe (e)=- 7 r- 1 In 1 [e-e fe -E(e)]- 1 . (19) 

In the paramagnetic state G is T-independent if we assume (S ■ er) = nS/2, and a depends on 
temperature only weakly through the Fermi function. If we neglect this thermal smearing around 
the Fermi energy p we may calculate at T = but consider the results to apply to the actual 
paramagnetic state at T > T c . We find [0] 

k 

where v k — Ve fc is the electron velocity and a 3 is the volume of the unit cell. Since A k depends 
on k only through e k we may define a function <fi(e) such that <j>'(e k ) = [^(m)] 2 - Hence the sum 
in equation ( |20| ) may be written as 

E Ve fc- v ^) = -E^) v2e fc' (2i) 

k k 

the last step following by means of Gauss's theorem. For a simple cubic tight-binding band e h = 
Y2r cos kpCi, with P summed over x,y,z, V 2 e fc = — a 2 e fc . Then the summand in equation ( |2l| ) 
is a function of e k only and equation ( |20| ) becomes 



2ne 2 f 

a = ^hJ dEEN c( E )^ E ) ( 22 ) 

where N (E) is the density of states for the simple cubic band. If N (E) is replaced by a suitably- 
scaled Gaussian N g (E) = (3/Tr) 1 / 2 ^-! 

exp [— 3(i?/VF) 2 ] , corresponding to an infinite dimensional 
approximation, equation (B2L) may be simplified by integrating by parts: 



ire 2 W 2 



dEN(E) [A E {,j)f . (23) 



Here A e (ji) is defined by the right-hand expression in equation (|T^) with e = /i, e k — E. 

Since it is convenient to use the elliptic density of states N (e) to calculate the Green function 
and self-energy, as in the previous section, it is reasonable to evaluate a using equation ( p2|) with 
-ZV C replaced by N Q . We take the strong-coupling limit J — > oo for simplicity and the results 
for the resistivity p = er -1 are plotted against band-filling n for various S in figure 2. We have 
taken a = 5 A which is comparable with the Mn — Mn spacing in manganites. For J = oo the 
band-width W is the only energy-scale and, since the integral in equation ( |22] ) is dimensionless, p 
does not depend on W. In fact in the DE regime JS > 2W the resistivity is almost independent 
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of both J and W. Use of the alternative formula (p^), with N replaced by N Q , gives very similar 
results. It is seen in figure ^| that p diverges correctly at n — 0, owing to the absence of carriers, 
and at n = 1 where the system becomes a Mott insulator. Well away from these insulating 
limits p does not depend strongly on S, so that quantum spin effects are not very important. 
Furthermore p w lmQ cm over a wide range of band-filling, which is much smaller than observed 
in some manganites above T c , as discussed in section 0. This agrees with the conclusion of Millis 
et al. that the DE model, with electrons scattered purely by disordered local spins, cannot 
describe the physics of the manganites completely. Early work by Furukawa Jl0| seemed to point 
to another conclusion, although the DMFT is equivalent to our CPA approach. We showed |j] 
that the confusion arose from Furukawa's use of a convenient, but rather unreasonable, Lorentzian 
density of states. Using a Lorentzian to replace N (E) and to calculate 4>(E) in equation (|22| ) 
leads to a divergent integral. However equation ( p3|) is similar to the form of a used by Furukawa 
and use of a Lorentzian to replace N c (E) and calculate A E {p) gives a convergent result. Results 
of such calculations for the limit J — > oo are shown in figure || and it is remarkable that p is at 
least an order of magnitude larger than one finds in figure || for the more reasonable elliptic band. 




0.5 - 



0.0 1 ' ' ' ' ' ' ' ' ' 1 

0.0 0.2 0.4 0.6 0.8 1.0 
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Figure 2: The zero field paramagnetic state resistivity p = a 1 versus band-filling n for the double- 
exchange model. Here J — oo, a = 5l, and S = 1/2, 1, 3/2, 2, 5/2 and oo, p increasing with S. 
The elliptical density of states and formula ( p2|) are used. 



5 Magnetism in the DE model 

As discussed at the beginning of section ||, the full equation of motion approach to many-body CPA 
is required to determine magnetic properties such as spin susceptibility \ an d Curie temperature 
T c . In our first paper this involved a hierarchy of Green functions satisfying AS + 1 coupled 
algebraic equations for local spin S; only the S = 1/2 case was briefly discussed. In a second paper 
a major simplification was achieved by introducing generating Green functions which generate 
all the required Green functions by differentiation with respect to a parameter. The coupled 
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Figure 3: As in figure but using a Lorentzian density of states and formula fl2 



equations are then replaced by a single first-order linear differential equation, the parameter being 
the independent variable, whose analytic solution yields the required CPA equations for the Green 
functions. The classical limit S = oo can then be taken and for J = oo the equations coincide 
with those of DMFT, which are only obtainable in the classical limit. Our many-body CPA is 
therefore an analytic approximation to DMFT for arbitrary quantum spin S which becomes exact 
for S = oo. The many-body CPA also coincides with Kubo's J|4| one-electron CPA in the limit 
n — ► where that is valid. 

To determine the magnetic properties one problem remains; the CPA and DMFT equations 
contain one set of correlation functions ((S z ) m ) which cannot be obtained directly from the Green 
functions. There is an indirect procedure for determining these correlation functions within CPA 
but it proves to be unsatisfactory, never yielding ferromagnetic solutions. However, for S — oo, 
DMFT provides a way to calculate the probability distribution function P(S Z ), and hence ((S z ) m ), 
and we used an empirical extension of this for finite S. This extension guarantees that the spin 
susceptibility exhibits the correct Curie laws for band occupations n = and n = 1. Thus for 
ii = we have a Curie law over the whole temperature range, corresponding to N independent 
spins S. For n = 1, with J = oo, we have independent spins S + ^- F° r < n < 1 wc find 
a finite Curie temperature T c and some results ^ are shown in figure ||. In figure 4(a) T c is 
plotted as a function of n for various S with J = oo, using the elliptic band. Clearly for finite 
S ferromagnetism is more stable for n > 0.5 than for n < 0.5, in agreement with the findings 
of Brunton and Edwards |2j| . For S = oo the result agrees closely with that of Furukawa |25| . 
In figure 4(b) we see the effect on T c for S — 1/2 of changing the bare band-shape from elliptic 
to simple cubic tight-binding. A dip in T c occurs around n — 0.3 which is the region where the 
ground state of the simple cubic DE model with S = 1/2 is rigorously not one of complete spin 
alignment J23|. 
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Figure 4: The Curie temperature k B T c /W of the double-exchange model versus band-filling n 
for various S, calculated with J = oo using the elliptic band (a); the effect on T c for S = 1/2 
of changing the elliptic band to the density of states for a simple cubic tight-binding band with 
nearest neighbour hopping (b). 

6 Many-body CPA for the Holstein-DE model 

This section is based on Green's recent study of the Holstein-DE model in which the electrons 
of the DE model couple to local phonons as in the Holstein treatment of small polarons |2(| p7| . 
The Hamiltonian is 



H = E vUv - J E s i ■ °i h E ( s i + <) 

ij a i i 

-sE n * ( & l+\)+"E 6 &- 



(24) 



The first three terms constitute the DE Hamiltonian of equation (|l|) while the first, fourth and 
fifth terms form the Holstein model. Einstein phonons on site i, with energy lo and creation 
operator b\, couple to the electron occupation number n i = n ia with coupling strength g. The 
electron-phonon coupling is of the form — g' J^. n i x i , where x i is the displacement of a shell of 
atoms surrounding site i, and in application to the manganites it may be regarded as an effective 
Jahn- Teller coupling. Previous studies of this model have either concentrated on coherent polaron 
bands, like Roder et al. pq] , or have treated the phonons classically |l^| so that there are no 
polaron bands at all. The many-body CPA approach is able to encompass both aspects and to 
describe the crossover from quantum polarons to the classical picture as temperature and/or model 
parameters are varied. The relationship to previous theoretical work and to experimental studies 
of the manganites is discussed fully in section ^. However we mention briefly below some related 
work on the pure Holstein model, without coupling to local spins. 

Sumi considered the Holstein model with one electron in the band, first treating the phonons 
classically [^9| and later quantum mechanically |50| . The classical case, with frozen displacements 
x v corresponds to a multicomponent alloy for which CPA is the best local approximation. In his 
dynamical CPA treatment of quantum phonons, Sumi p0| treated the one-site dynamics correctly 
and his work is completely equivalent to the more recent DMFT treatment of Ciuchi et al. [ pl| . 
As a general rule dynamical CPA and DMFT are the same for one-electron problems. DMFT 
is the correct extension of CPA to the many-body problem of finite electron density but for the 



11 



Holstcin model, as for the DE model, it cannot be carried through analytically in the quantum 
case. Numerical work [[32], applying DMFT to the Holstein model has been aimed mostly 
at understanding superconducting transition temperatures and charge-density-wave instabilities 
rather than the polaron physics we are mainly concerned with. An unfortunate feature of the 
Holstein model for spin 1/2 electrons is that in a quantum treatment the true ground state for 
strong electron-phonon coupling consists of unphysical singlet bipolarons with two electrons bound 
on the same site. This problem does not occur in the one-band Holstein-DE model since strong 
coupling J to local spins prevents double occupation of sites, as pointed out earlier. It is also 
bypassed if the phonons are treated classically, as in the work of Millis et al. J|j| on the Holstein 
model. The Holstein model is more complicated than the DE model and it turns out that our 
many-body CPA no longer reduces to the correct one-electron dynamical CPA/DMFT |3(], |3lJ as 
band-filling n — > 0. Although correct in the atomic limit = 0, our theory is clearly cruder for 
the Holstein and Holstein-DE model than for the pure DE model. 

We start by deriving the Green function for the Holstein-DE model in the atomic limit. The 
Hamiltonian H AL in this limit is given by equation (^ij) with the first term omitted and with 
site indices and summation suppressed. We remove the electron-phonon coupling by the standard 
canonical transformation Q] H = e s H AL e~ s where s — —(g/u>)n(b^ — b). Under this transforma- 
tion b — > b + (g/ui)n and the Hamiltonian separates into a fermionic and bosonic component: 

H = H f + H h (25) 



H { = -JS -<r-h{S z + a z )~ (g 2 /uj) n 2 , H h = uob ] b . (26) 

The transformation corresponds to a displacement of the equilibrium position of the phonon 
harmonic oscillator in the presence of an electron and the downward energy shift g 2 /w is a polaron 
binding energy which we write as Aw, where A = g 2 jui 2 . If two electrons occupy the site (n = 2), 
which will not occur for large J, the energy shift becomes 4g 2 /ui 2 corresponding to an on-site 
bipolaron. Writing out explicitly the thermal average in the definition of the one-particle retarded 
Green function we have 



Gf(i) = -i0(*)(k(t),4] + ) 



Tr{e-^A L [ C(T (t), c t] + } 



= ^ Trie^} (27) 

and the canonical transformation introduced above can be carried out within the traces, using the 
property of cyclic invariance. Thus H AL — > H , c£ — > ^ c l an d c (i) becomes 

e im Xc„e~ im (28) 

where X — exp \g(tf — b)/uf\. Using equation (p5|), we can write the traces in equation j27| ) as 
products of fermionic and bosonic traces. Hence we find 

Gf (t) = -i9(t) {<c CT (t)ct) f F(t) + <4 C(7 (t)> f F*(t)} (29) 

where F(t) = (X(t)X^) b and the thermal averages ( )t,( ) b correspond to the systems with 
Hamiltonians H i and H h respectively. It may be shown j34| that 

F(t) = c - x(2b+1 '> cxp {2A [b(b + 1)] 1/2 cos [u(t + i/9/2)]} (30) 

where b = b(u>) = (c 13 ^ — l) 1 is the Bose function with (3 — {k B T) 1 . The last factor is of the 
form cxp (zcos</>) which generates the modified Bessel functions I r (z): 

oo 

exp (z cos 0) - U z )^- (31) 
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To evaluate the fermionic averages we consider for simplicity the limit J — > oo in zero field (h = 0). 
Then the last term in H f may be written — (g 2 /uj)n, since n — or 1 only, and this may be absorbed 
into the chemical potential which is finally determined to give the correct number of electrons n per 
atom. Thus H { is just the DE Hamiltonian in the atomic limit and the sum of the two fermionic 
averages corresponds to the function G AL (t) whose Fourier transform is given by equation (|l6|). 
It is easy to see that the first and second thermal averages in equation ( |29| ) take constant values 
(1 - n)(S +1)/(2S +1) and n(S + l/2)(25+ 1) respectively. Hence, from equations @-(|l|), we 
obtain the Fourier transform of G AL , with J — > oo and h = 0, in the form 

r AL M = y U2A [&M(&M + 1)] 1/2 } (gg + l)|e^ + (S + 1)(1 - n)e~^ 

{£> r f^ o (25+l)exp{A[26H + l]} e + ru; ' 1 J 

The density of states — 7r _1 Im G AL (e) is shown in figure || for the classical spin limit S — ► oo 




Figure 5: One-electron spectra of the Holstcin-DE model in the atomic limit at zero and very high 
temperature. They consist of delta-functions, with energy spacing to, whose strength is indicated 
by the envelope curves. The plots are for the paramagnetic state with S = J = oo, h = 0, 
n = 0.5, lo/W = 0.05 and g/W — 0.18, where W is a unit of energy later to be identified with the 
half-width of the electron band in the full Hamiltonian. 

at quarter-filling n = 0.5. It consists of delta-function peaks separated in energy by ui and the 
envelope curves show the weight distribution at low and high temperature. W is an energy unit 
which, when we go beyond the atomic limit, will be the half-width of the itinerant electron band, 
as usual. The values adopted for the parameters w/W and g/W relate to the manganites, as 
discussed in section 0. The symmetry of the spectrum about zero energy is due to the choice 
of filling n = 0.5; in general at T = the lower and upper "bands" have weights n and 1 — n 
respectively. By counting weights it may be seen that for any n the chemical potential lies in the 
peak at e = 0, which has very small weight e _A /2 per spin. The shape of the envelope function 
at T = 0, with two maxima and very small values at the centre of the pseudogap between them, 
may be understood physically as follows. The delta-function at e = rw (r > 0) corresponds 
to an excitation from the ground state, with no electron and the undisplaced oscillator in its 
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ground state, to a state with one electron and the displaced oscillator in its r th excited state. The 
strength of the delta- function depends on the square of the overlap integral between the displaced 
and undisplaced oscillator wave functions. Clearly this is very small for r = and goes through 
a maximum with increasing r as the normalized displaced wave function spreads out. At T = 



it is easily seen from equation (32), using l r (z) ~ (z/2)' r l/|r|! for small z, that the weight of the 
delta-functions at e = ±rw is proportional to A r /r! . Hence the maxima in the envelope curve 
occur at e w ±Xlo, which is the polaron binding energy. 

We now turn to the Holstein-DE model with finite band-width. As for the DE model it is 
necessary to use the full equation of motion method to derive the many-body CPA in the presence 
of a magnetic field and/or magnetic order j|. In the present case it is very difficult to determine 
self-consistently all the expectation values which appear. We therefore approximate them by their 
values in the atomic limit. It then turns out that in the zero field paramagnetic state, for J = oo 
and with the elliptic band, the CPA Green function G again satisfies equation (|l4[), with G AL now 
given by equation ([32]). 

The densities of states calculated for T — using equations (|32] ) and Jl4| ) with S = oo, n = 0.5, 
u>/W — 0.05 and various values of g/W are shown in figure q, Apart from lacking the perfect 
symmetry about the chemical potential /i — the results are qualitatively similar for other values 
of n not too close to or 1. For g = we recover the elliptic band with half- width W/y/2 as for the 
DE model with J = oo, S — oo. As g increases the density of states broadens and small subbands 
are split off from the band edges. As g increases further a pseudogap develops near the chemical 
potential. At a critical value g — g c a gap appears which contains a small polaron band around the 
chemical potential. Increasing g further causes more bands to be formed in the gap, with weights 
similar to those of the relevant atomic limit. It should be pointed out that the paramagnetic 
state considered here at T = is not the actual ground state, which is ferromagnetic. We discuss 
the magnetic state later. The effect of increasing temperature on the density of states in the gap 
region is shown in figure [7] for g = 0.18W > g c . With increasing T the polaron bands grow rapidly 
and eventually merge to fill the gap. 

It is important to compare these results with the standard small polaron theory developed by 
Holstein [25, 33]. Holstein distinguished between "diagonal transitions", in which the number of 
phonons is unchanged as the electron moves from site to site, and "nondiagonal transitions" in 
which phonon occupation numbers change. The former give rise to a coherent Bloch-like polaron 
band of half-width We~ x ^ 2b+1 ^ which decreases with increasing temperature. The nondiagonal 
transitions are inelastic processes which destroy phase coherence and the polaron moves by diffusive 
hopping. The hopping probability increases with temperature so that polaron motion crosses over 
from coherent Bloch-like at T — to diffusive hopping as k B T approaches the phonon energy 
u>. The paramagnetic state of the Holstein-DE model differs from this standard picture in one 
important respect. There are no well-defined Bloch states, owing to strong scattering by the 
disordered local spins, so no coherent polaron band will form. This is fortunate because the CPA 
treatment of electron-phonon scattering will never lead to coherent states of infinite lifetime at the 
Fermi surface at T — 0. However in the presence of strong spin disorder it should be satisfactory. 
We interpret the central band around the chemical potential in figure |?] as an incoherent polaron 
band whose increasing width as the temperature rises is due to life-time broadening of the atomic 
level. The life-time decreases as the hopping probability increases with rising temperature. 

To substantiate this picture we study the central polaron band in the limit of very strong 
electron-phonon coupling. In this limit it can be shown that we need retain only the r = term 
in equation (32) and it is then easy to solve equation dl4| ) for G. The result is of the same form 
as equation (17) but with 

D * = J. W 2 e -H2b+l) lo [b{b + 1}] l/2^ (33) 

The central band is thus elliptical with half- width D and weight D 2 /W 2 . It is now easy to calculate 
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Figure 6: The one-electron density of states (DOS) for the Holstein-DE model with half-bandwidth 
W, for the hypothetical paramagnetic state at T = 0, with various strengths of electron-phonon 
coupling g/W. Other parameters as in figure ||. 




Figure 7: Evolution with temperature f3 — (k B T) of the polaron subbands in the pseudogap 
around the chemical potential /i = for g/W = 0.18. These subbands at T — can just be seen 
in figure 0. All parameters as in figure ||. 
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the conductivity a from equation ( p0| ) and, using D 2 <C W 2 , we find 
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The last step follows by using the asymptotic forms for strong coupling and high temperature 
I (z) ~ (2irz) 1 ^ 2 expz and b ~ 1 . The temperature dependence of er is the same as in 

the standard theory of small polaron hopping conduction but with activation energy Aw/4 
equal to one quarter, instead of one half, of the polaron binding energy. The standard result is 
for one electron coupled to the lattice without spin disorder. Nevertheless this establishes the 
link between our work and standard small polaron theory in the strong coupling limit. However 
the results shown in figure Q with parameters relevant to typical manganites, are far from this 
limit. They correspond to intermediate coupling and in the actual paramagnetic state above the 
Curie temperature the polaron bands are largely washed out. In this regime, with increasing 
temperature, there is a crossover from polaronic behaviour to a situation where the phonons 
behave classically, the case considered by Millis et al. Q. For electron-phonon coupling greater 
than a critical value these authors find a gap in the density of states which gradually fills with 
increasing temperature. However in their classical treatment there are no polaron bands in the 
gap so that the link with standard polaron physics is not established. 

Apart from the symmetry about e = the above results for n — 0.5 are not untypical of the 
general case. For general n the main lower and upper bands, separated by a gap for g > g c , have 
approximate weights n and 1 — n respectively. The chemical potential at T — is always confined 
to the polaron band arising from the r = term of equation (|3^), and moves from the bottom at 
n = to the top at n = 1, so that we correctly have an insulator in these limits. 

To calculate the Curie temperature T c we need the full CPA theory combined with an exact 
result of DMFT for S = oo |J . Results on T c for the same parameters as before are plotted as 
functions of electron-phonon coupling g in figure g[ The suppression of T c with increasing g was 



first noted by Roder et al. |28 and our own results are quite similar to those of Millis et al. [[L2| 
In our CPA we have no reliable means of calculating the probability distribution function P(S Z ), 
so to go below T c we use the mean-field approximation for the ferromagnetic Heisenberg model 
with classical spins and nearest neighbour exchange. The exchange constant is determined by T c . 
We plot the up- and down-spin density of states for T = 0.005W/fc B <C T c and g = 0.16W > g c 
in figure ^, also showing curves for the saturated ferromagnetic state and paramagnetic state at 
T = for comparison. The value g = 0.161^ is closer to g c than the value of 0.18W used in 
figures and ^ and we discuss these results in relation to the manganites in the next section. In 
figure Owe plot the resistivity p as a function of temperature, for the same parameter set, with 
different applied fields h. The resistivity peaks sharply at T c , a nd for comparison we show results 
for weaker electron-phonon coupling g/W — 0.10 in figure P. The curve in figure [ll] is almost 
indistinguishable from that of figure 7 in reference || for g/W = 0.01. This is not surprising 
since we see from figures |6| and ^ of this paper that the density of states and T c change very little 
between g/W — and g/W = 0.1. These results are all discussed further in the next section. 



7 Application to manganites 

Much experimental work has concentrated on the systems La 1 _ x Ca x Mn0 3 and La 1 _ x Sr x Mn0 3 
with x « 0.33, in the middle of the ferromagnetic regime where CMR is observed. For brevity 
we denote these systems by LCMO and LSMO respectively. In comparing experiments with the 
detailed results of section || it should be borne in mind that the deduced parameters will be 
influenced by our convenient choice of n = 1 — x = 0.5 rather than n = 0.6 — 0.7. However the 
correct general picture should emerge. Since we consider a homogeneous state we are not concerned 
with the existence of charge ordering for x = 0.5. We should also discuss the effect of using a 
one-band model for the e g band, rather than the more realistic two-band model. An important 
point to notice is that in the two-band model large Hund's rule coupling J to the local spins 
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Figure 8: Suppression of the Curie temperature of the Holstein-DE model with increasing electron- 
phonon coupling g/W. The plot is for S = J = oo, h = 0, n = 0.5 and w/W — 0.05. 
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Figure 9: The up- and down-spin density of states of the Holstein-DE model with g/W — 0.16 
for k B T = 0.005FF < k B T c where (S z ) = 0.915. Also shown are the DOS for the saturated 
ferromagnetic state at T = and for the hypothetical paramagnetic state at T = 0. All plots are 
for S = J = oo, h = 0, n = 0.5, lu/W = 0.05 and g/W = 0.16. 
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Figure 10: Resistivity p versus temperature for the Holstein-DE model with S = J = oo, n = 0.5, 
uj/W = 0.05, intermediate coupling g/W = 0.16 and various applied fields h. The lattice constant 
is taken as a = 5A, slightly larger than the Mn — Mn spacing in the manganites. 
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Figure 11: The same plot as figure but for weak electron-phonon coupling g/W = 0.10. 
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no longer produces a Mott insulator for n — 1 and it is necessary to introduce on-site Coulomb 
interaction |55|, p6[ . Millis et al. |l2| do not do this so that for n = 1, which should correspond to 
the undoped insulator, they have a metal. Also, for weak electron-phonon interaction, they find 
T c is largest for this value of n. This contrasts strongly with figure ^ where, in the one-band DE 
model, T c = at n — and 1 and has a maximum in between. Held and Vollhardt p6| , using 
DMFT, obtain very similar results to ours for S — oo in the two-band model with strong on-site 
Coulomb interaction included. However their values of T c are about twice ours. Nevertheless 
we conclude that the one-band DE and Holstein-DE models may sensibly be used to discuss the 
manganites. We see from figures || and ^ that for the DE model neither T c nor the resistivity p 
vary enormously with S so that S — oo is a reasonable approximation to the 5 = 3/2 Mn spin. 
It is also reasonable to take the DE limit J — > oo |37|j . For comparison with experiment we take 
W = 1 eV which band calculations ]|$| [5?j suggest as an appropriate half-bandwidth for the e g 
band. Also in figures 6-11 we have taken lo/W = 0.05 to correspond to observed transverse optic 
phonons with to sa 40 — 70meV which couple strongly to the electrons in LCMO |59) . 

Perhaps the most striking feature of the manganites is the very different behaviour observed in 
apparently similar materials such as LSMO and LCMO. For LSMO, with x w 0.33, T c « 370 K 
whereas for LCMO, with a similar x, T c « 240 K. The difference in behaviour of the resistivity p 
above T c is much more striking. For LSMO p«4 mfl cm and increases slowly with temperature 
as in a poor metal |4(J. The p{T) curve is very similar to that of figure [y] for g/W — 0.1 except 
for a much larger resistivity at low temperature in our calculations. Since this feature persists 
even for g/W = 0.01 (figure 7 in reference |^]) it presumably arises from overestimated spin 
disorder scattering at low temperatures due to our use of the classical spin Heisenberg model to 
determine P(S Z ). In LCMO the resistivity rises to a maximum at T c of about 40 mfl cm and then 
falls with increasing temperature above T c [JO]. In contrast to LSMO there is thus a transition 
from metallic to insulating behaviour. Also the resistivity peak is strongly reduced and shifted to 
higher temperature with increasing applied magnetic field. This is the CMR effect. This type of 
behaviour is seen in figure for g/W — 0.16. The main differences between theory and experiment 
are a more rapid observed drop in p with decreasing temperature below T c and a more sensitive 
observed CMR effect. h/W = 0.004 corresponds to a field of about 20 T for W = 1 eV and the 
corresponding reduction in p in figure [l^ is achieved with a field of about 5 T experimentally. 
Millis et al. || noted a similar problem in their work using classical phonons. Both of the 
discrepancies mentioned might be remedied by introducing a dependence of g on p, corresponding 
to more efficient screening of the electron-phonon interaction with increasing metallization. The 
huge reduction in resistivity peak on reducing from g/W from 0.16 to 0.10 shows the extreme 
sensitivity of p to changes in g. The main point to notice is that we can understand the enormous 
difference between LCMO and LSMO within our theory by assuming the electron-phonon coupling 
changes from g/W = 0.16 in LCMO to g/W = 0.10, or slightly greater, in LSMO. The observed 
ratio of the Curie temperatures, slightly less than 2, is then in accord with figure ||. As discussed 
in section || the critical coupling g c for the formation of a polaron band is g c /W « 0.15, with our 
phonon energy lo/W = 0.05, and to obtain the right order of magnitude for p above T c in LCMO, 
g is pinned down closely to 0.16. A larger value for g leads to too high a resistivity and too low a 
Curie temperature. It is interesting that neither p(T) nor T c change when g/W is varied between 
0.1 and 0. This means that LSMO can be described very well by the pure DE model, as stressed 
by Furukawa [ p5| , but it does not follow that electron-phonon coupling is negligible. However, 
from the results of Millis et al. jl2| for classical phonons, one can understand why a coupling small 
enough to give a LSMO-like p{T) curve does not lead to a change in slope of the rms oxygen 
displacements, as a function of temperature, at T c . No such change is found in LSMO [Q, in 
contrast to the case of LCMO . It is more difficult to understand the observation [Q of static 
local Jahn- Teller distortions in LSMO at room temperature, apparently associated with localized 
carriers in the presence of metallic conduction. 

From figures || and || we see that for g/W = 0.1, appropriate to LSMO, there is no sign of a 
pseudogap in the density of states. An actual gap in the hypothetical paramagnetic state at T = 
appears at g — g c with g c between 0.15 and 0.16. From figure ^| we see that for g/W — 0.16, 
appropriate to LCMO, a few polaron subbands have appeared in the gap. These are seen much 
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more clearly in figure f^for g/W — 0.18 when there is a larger gap. However the subband structure 
is washed out completely for (3W = 25, corresponding to T = 464 K for W = leV, and this effect 
will occur at a much lower temperature for g/W = 0.16. Thus in the actual paramagnetic state 
of LCMO above T c we do not expect the quantum nature of phonons to manifest itself, so we 
are essentially in the classical regime of Millis et al. p2] |. The same is true in the saturated 
ferromagnetic state at T — where only a pseudogap appears in figure ^. As the temperature 
rises towards T c a minority spin band grows, also with a pseudogap, while the majority spin 
band loses weight. The width of the ferromagnetic bands decreases with increasing temperature, 
corresponding to the DE effect, but the narrowing in the paramagnetic state is not so marked 
as in the pure DE model. Thus double exchange is not so effective in the presence of strong 
electron-phonon coupling, which is consistent with the reduction in T c shown in figure ||. 

Clearly our picture of the manganites is close in spirit to that of Millis et al. 12 , although the 
relationship to polaron physics is not so clear in their classical approximation. Other authors adopt 
completely different viewpoints. Furukawa p5f rejects the importance of electron-phonon coupling 
and argues in favour of phase separation models of LCMO with, for example, ferromagnetic and 
charge-ordered regions. Such scenarios have been extensively discussed by the Florida group 
|45[ f46f . It is controversial whether such phase separation can occur for x w 0.33 which is far 
from the antiferromagnetic insulating phase near x = and the region with charge and orbital 
ordering near x = 0.5. Nagaev J47I also argues against any polaronic effects, while Alexandrov 



and Bratkowsky 1 48 , [19| assume strong electron-phonon coupling with small polarons even in the 
ferromagnetic state and with immobile bipolarons forming near T c . The public correspondence [50| , 
fH^ I between Nagaev and Alexandrov and Bratkowsky (AB) centres on estimating the magnitude 
of the polaron binding energy E and the criterion for small polaron formation | pl| . Since our 
picture of LCMO lies between their extreme views it is interesting to compare our estimates with 
theirs. For LCMO we find E = g 2 /oj ~ 0.5 eV for W — 1 eV whereas, for manganites in general, 
Nagaev estimates £ p k 0.1 - 0.3 eV and AB estimate E ~ 1 eV. Our condition for small-polaron 
formation in a paramagnetic state at T = is g > g c « 0.15VF which corresponds to E p > 0A5W. 
Nagaev adopts the criterion E p > W, remembering that W is the half-bandwidth in our notation, 
while AB |5l| propose E p > 2W(8z)~ 1 / 2 = 0.29W with number of nearest neighbours z taken 
as 6. AB's condition is less stringent than Eagle's condition for "nearly small polarons" 
E > Wz- 1 / 2 = 0A1W which is close to ours. Both as regards this criterion and the value 
of E for LCMO, our results are intermediate between Nagaev's and AB's, as expected. For 
LSMO, on the other hand, our estimate of E p is 0.2 eV. In this case we agree with Nagaev that 
electron-phonon coupling is not important. AB assume stronger electron-phonon coupling than 
we find in LCMO. Hence they have a coherent small-polaron band in the ferromagnetic state 
which apparently supplies the required metallic conductivity. The carriers are unable to bind to 
form singlet bipolarons. 

AB [[l8| [l9| assume that the strong Hund's rule coupling, which we and most authors assume 
is effectively infinite, can be treated in mean field theory. Instead of the double-exchange effect, 
in which the occupied polaron band merely narrows as its weight redistributes between up and 
down spin, the Hund's rule exchange splitting collapses above T c and it is assumed that immobile 
singlet bipolarons can form. Hence near T c there is an enormous increase in resistivity which 
decreases above T c as the bipolarons dissociate. The problem with this theory as it stands, even 
more serious than doubts about the existence of polarons above and below T c , is the mean field 
treatment of the large Hund's rule exchange. 

The giant isotope effect observed in LCMO and Nd 7 Sr 3 Mn0 3 has been interpreted as ev- 
idence for immobile bipolaron formation |p3[ . The effect is seen as a decrease of T c and a big 
increase in resistivity, particularly near T c , when 16 O is replaced by 18 O. To investigate this effect 
in our model some care is needed. It was pointed out that in equation ( p4] ) the electron-phonon 
coupling term corresponds to a term —g' n i x i where x i is to be associated with oxygen dis- 
placement around a Mn atom. Here g' should be independent of the oxygen mass M. In the 
second-quantized form of equation ( pi]) one finds g — g' (2Mlo) 1 ^ 2 so that the polaron binding 
energy g 2 /ui — g' / (2Moj 2 ). Since for an oscillator oj cx M^ 1 / 2 the polaron binding energy g 2 / 
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Figure 12: The absence of an isotope effect in LCMO due to a simple mass scaling (oc M 1 / 2 ) of 
the phonon frequency. 



is expected to be independent of M. However g varies as Af^ 1 / 4 . We have therefore recalculated 
the resistivity p{T) with the same parameters as used in figure^ for 16 but with g and u> scaled 
appropriately for 18 O. The results are compared in figure The almost complete absence of an 
isotope effect in both T c and p is remarkable in view of the moderately strong electron-phonon 
coupling in our model of LCMO. Since we do not believe in the bipolaron theory for LCMO 
another explanation of the isotope effect must be sought. Nagaev has considered several possi- 
bilities, all based on the isotope dependence of the number of excess or deficient oxygen atoms 
in thermodynamic equilibrium p7| , p4| . In LCMO with x — 0.2 there is experimental evidence 
that at least part of the isotope effect is due to varying oxygen content J55| . A change in oxygen 
nonstoichiometry could lead to a change in carrier density and, perhaps due to volume change, 
to a change in electron-phonon coupling g' . In view of the sensitivity of p(T) to electron-phonon 
coupling we think that the isotope effect may be due more to a change in g' , or to a change in 
uj beyond the simple mass scaling, than to a change in n. Since Nagaev believes electron-phonon 
coupling to be unimportant he did not consider this particular consequence of his nonstoichiometry 
proposal. The involvement of electron-phonon coupling in the isotope effect would explain why 
the Tq shift in LSMO is much smaller than in LCMO |^6| . To illustrate such a mechanism we have 
recalculated p(T) for the 18 O system with scaled u> as before (oc AT -1 / 2 ) but with g unchanged 
from its 16 value. This implies a 3% increase in g' . The results are shown in figure 13 
the paramagnetic state they are quite similar to those observed in LCMO with x 
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Nd Q 7 Sr Q 3 Mn Q 3. The only conclusion we can draw from this arbitrary calculation is that it may 
be extremely difficult to determine the true origin of the isotope effect in such a critical system 
as LCMO on the verge of small-polaron formation. However it does not arise from a simple mass 
scaling of the phonon frequency. 

Another example of the sensitivity of p to variation in parameters is the effect of pressure [ pT| . 
The strong suppression of the resistivity peak and the increase in Curie temperature in LCMO 
can be modelled within our theory by increasing the bandwidth [^|, and assuming other terms 
in the Hamiltonian are constant. Calculated results are shown in figure This pressure effect 
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Figure 13: A giant isotope effect due to mass scaling of the phonon frequency together with a 3% 
increase of electron-phonon coupling strength in the 18 O system. This arbitrary modelling shows 
how an effect similar to that observed ]53| might arise from small parameter changes associated 
with changed oxygen content. 

is due to a reduction in g 2 /(ujW), the ratio of polaron binding energy to the half-bandwidth. 
A simple estimate, using the known compressibility and dependence of W on lattice constant, 
shows that the theoretical pressure for a given effect is about four times larger than that required 
experimentally. This is the same factor that we found in the case of the magnetic field required for 
a given CMR effect, so both discrepancies could possibly be removed with the same dependence 
of g on p , due to screening, postulated earlier. 

According to our theory we expect pseudogaps at the Fermi level to be observable in the 
density of states of LCMO both below and above T c . These should appear in experiments such 
as scanning tunnelling spectroscopy, photoemission and optical conductivity measurements. No 
pseudo gaps are expected in LSMO. The pseudogap is a feature of the atomic limit, typified 
by the envelope function in figure ^| for T = with maxima determined by the polaron binding 
energy g 2 /w, and is completely washed out when g 2 /(u>W) < 0.2 (see figure ||). Early results 
of scanning tunnelling spectroscopy on LCMO with x = 0.3 seem unlikely to relate to the 
bulk. In the ferromagnetic state at 77 K there is a huge gap of about 1 eV. ft is not clear why the 
authors interpret this as evidence for half-metallic ferromagnetism. A gap of this size associated 
with small-polaron formation in the bulk would imply an unrealistically large electron-phonon 
coupling, certainly incompatible with metallic conduction and a Curie temperature of reasonable 
magnitude for LCMO. More recently Biswas et al. [^| reported a scanning tunnelling spectroscopy 
study of several manganites. The results are very much in accord with our theory. There is no 
gap in the low temperature ferromagnetic state but a small gap (pseudogap) appears for T « T c 
in the low (high) T c materials. As T increases above T c the pseudogap or gap gets filled in as we 
would expect (see figure 0) . 

In an extremely interesting paper on angle-resolved photoemission spectroscopy for the bilayer 
manganite La-L 2 Sr 1 g Mn 2 7 , nominally with n = 0.6, Dessau et al. [ p0| interpret their results very 
much in the spirit of our theory. In this layered structure we cannot estimate the strength of 
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Figure 14: The effect of pressure on the resistivity and Curie temperature in LCMO by increasing 
the half-bandwidth W to W = 1.05W and 1.1W. 



the electron-phonon coupling from T c , whose low value of 126 K is largely a result of quasi- two 
dimensional fluctuations. However the large low-temperature resistivity of more than 3mf2cm 
suggests that small-polaron bands might exist in both the ferromagnetic and paramagnetic state. 
This indicates a larger value of g/W than in the cubic manganites, which might be expected 
from a reduced bandwidth in the layered structure. In this case the polaron bands would not be 
washed out above T c and should exist within a wider gap than in the ferromagnetic state. This 
is compatible with the observations. Dessau et al. |30| interpret the widths of their fe-resolved 
spectral peaks in terms of multiple polaron subbands. The fact that strong peaks all occur at 
0.65 — 1.0 eV below the Fermi level, independent of fc, suggests that the system is not so far from 
the atomic limit, with the spectral peaks given by the envelope function in figure for T = 0. 
This indicates a large polaron binding energy of 0.65 — 1.0 eV and Alexandrov and Bratkovsky 
have applied their theory [|TJ to these bilayer manganites. They interpret the observed upward 
shift in frequency v of the maximum in the optical conductivity <j{v), between 10 K and just above 
T c , as due to bipolaron formation. We cannot discuss the possible binding of two polarons on 
nearest-neighbour sites within our local approximation. However one docs not need bipolarons to 
understand an upward shift of the optical conductivity maximum with increasing temperature. 
This is clearly seen already in the calculations of Millis et al. Jl^] where the phonons are treated 
classically. In connection with the work of Dessau et al. , it should be mentioned that Moreo 
et al. |]6lf interpret the observed pseudogap not as an intrinsic property but in terms of phase 
separation. 

Kim et al. |^2| have reported very interesting measurements of optical conductivity a(y) in 
La 7 _ y Pr y Ca 3 Mn0 3 for y = 0, 0.13, 0.4 and 0.5. The Curie temperature takes corresponding 
values of approximately 245, 240, 155 and 120 K. With our parameterization this corresponds, 
from figure ||, to g/W varying between 0.16 and 0.22. The observed variation of o{v) with y and 
temperature T, over a wide range of photon energy up to 2 eV, is in general agreement with the 
classical treatment of Millis et al. |l^| . It is possible that the unusual behaviour in the far infrared 
region below 0.15eV, for y = 0.4 and 0.5, may relate to our calculated electronic structure in 
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the pseudogap in addition to direct excitation of phonons. Detailed calculations of angle-resolved 
photoemission and optical conductivity spectra, and their comparison with experimental data, 
will be reported later. 

8 Conclusion 

We have summarized the many-body CPA treatment of the double-exchange model, with and 
without coupling to local phonons, which is presented in detail in three earlier papers [§, ^ Q]- 
The method can deal with quantum local spins and phonons and is equivalent to dynamical mean 
field theory in the classical limit. For quite a wide range of fairly weak electron-phonon coupling 
we find that the results are essentially those of the pure double- exchange model. As the electron- 
phonon coupling is increased our theory can describe, for the first time, the crossover from the 
classical phonon limit fl2]| to the formation of small-polaron bands, and finally it links up with 
standard small-polaron theory in the strong-coupling limit. In section ^ we have given a much fuller 
discussion of the application to manganites than in the earlier papers. As well as basic properties 
such as Curie temperature, resistivity and magnetoresistance, we discuss the interpretation of the 
giant isotope effect, pressure effects and pseudogaps observed in scanning tunnelling spectroscopy 
photoemission and optical conductivity. 

We find that a typical manganite like La 1 _ x Ca x Mn0 3 , with x w 0.33, is in the critical regime 
on the verge of small-polaron formation. This explains its extreme sensitivity to changes in 
parameters arising from isotopic substitution and pressure. On the other hand La 1 _ x Sr x Mn0 3 , 
with weaker electron-phonon coupling, is quite well described by the pure double-exchange model, 
as pointed out by Furukawa |2j| . The large pseudogap seen in a bilayer manganite |30| indicates 
stronger electron-phonon coupling than in most cubic manganites. 

We are grateful to Lesley Cohen, D.M. Eagles, K. Kamenev, N. Furukawa and K. Kubo for 
helpful discussion and to the UK Engineering and Physical Sciences Research Council (EPSRC) 
for financial support. Elizabeth Rowley has given valuable help in preparing the manuscript for 
publication. 
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